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Improving  imaging  through  turbulence  via  aperture  partitioning 


Brandoch  Calef 

Boeing  LTS  Maui ,  535  Lipoa  Parkway  Suite  200,  Kihei  HI  96753 

ABSTRACT 

Speckle  imaging  techniques  make  it  possible  to  do  high-resolution  imaging  through  the  turbulent  atmosphere  by  collecting 
and  processing  a  large  number  of  short-exposure  frames,  each  of  which  effectively  freezes  the  atmosphere.  In  severe  seeing 
conditions,  when  the  characteristic  scale  of  atmospheric  fluctuations  is  much  smaller  than  the  diameter  of  the  telescope, 
the  reconstructed  image  is  dominated  by  “turbulence  noise”  caused  by  redundant  baselines  in  the  pupil.  I  describe  a 
generalization  of  aperture  masking  interferometry  that  dramatically  improves  imaging  performance  in  this  regime.  The 
approach  is  to  partition  the  aperture  into  annuli,  form  the  bispectra  of  the  focal  plane  images  formed  from  each  annulus, 
and  recombine  them  into  a  synthesized  bispectrum  from  which  the  object  may  be  retrieved.  This  may  be  implemented 
using  multiple  cameras  and  special  mirrors,  or  with  a  single  camera  and  a  suitable  pupil  phase  mask.  I  report  results  from 
simulations  as  well  as  experimental  results  using  telescopes  at  the  Air  Force  Research  Lab’s  Maui  Space  Surveillance  Site. 

1.  INTRODUCTION 

When  images  are  collected  through  the  atmosphere,  they  are  distorted  due  to  refractive  index  fluctuations  along  the  tur¬ 
bulent  path  [1].  The  object  may  be  recovered  by  the  bispectrum  technique,  which  exploits  phase  closure  and  atmospheric 
statistics  to  provide  an  estimate  of  the  object’s  Fourier  phase.  As  the  ratio  D/ro  of  the  telescope  aperture  diameter  to 
the  atmospheric  seeing  cell  size  increases,  the  optical  transfer  function  (OTF)  is  attenuated  and  the  number  of  redundant 
baselines  in  the  pupil  increases,  driving  down  the  SNR  of  the  bispectrum  estimate. 

The  usual  way  of  overcoming  the  latter  degradation  (which  is  sometimes  called  “baseline  redundancy  noise”)  is  to 
collect  short  exposure  images  for  many  atmospheric  decorrelation  times.  Averaging  together  their  bispectra  increases 
the  SNR  of  the  object  bispectrum  estimate  in  proportion  to  the  square  root  of  the  number  of  independent  atmospheric 
realizations.  However,  this  approach  is  not  always  practical.  For  example,  suppose  that  the  subject  of  observation  is 
a  satellite  passing  overhead.  Every  few  seconds,  its  orientation  relative  to  the  observer  changes  by  an  amount  that  is 
significant  at  the  resolution  of  the  telescope.  This  puts  a  limit  on  the  time  interval  of  frames  that  can  be  safely  co-processed. 

Another  way  of  defeating  redundancy  noise  is  to  mask  the  aperture,  either  with  a  pattern  of  holes  whose  diameters  are 
on  the  order  of  ro  (a  non-redundant  mask)  or  with  a  thin  annulus  at  the  outer  edge  of  the  pupil  (a  partially  redundant  mask). 
This  approach  is  called  aperture  masking  interferometry  [2] .  It  eliminates  or  greatly  reduces  baseline  redundancy  noise  at 
the  expense  of  a  much  lower  signal.  This  means  that  it  can  only  be  applied  when  the  object  is  bright  or  when  frames  can 
be  collected  over  a  long  period  of  time  to  accumulate  signal. 

In  the  following  pages,  I  describe  a  way  to  extend  this  technique  by  partitioning  the  pupil  and  measuring  the  speckle 
patterns  associated  with  each  section.  In  this  way,  the  redundancy  noise  is  reduced  while  still  making  use  of  all  light 
incident  on  the  aperture.  The  natural  partition  is  a  set  of  concentric  annuli,  with  the  outermost  ring  providing  information 
up  to  the  diffraction  limit  of  the  telescope. 

2.  BISPECTRUM  RECONSTRUCTION 


2.1  Review  of  bispectrum  technique 

The  object  will  be  reconstructed  from  speckle  images  using  the  bispectrum,  which  we  now  review.  Let  i  (v)  be  a  speckle 
image  (focal-plane  intensity  pattern)  formed  by  convolving  an  object  o{x)  with  an  atmospheric  point- spread  function 
(PSF)  h(x).  In  Fourier  space,  this  is  written  I(u)  =  0(u)H(u).  The  bispectrum  of  i(x)  is  defined  by  Bi(u,  v)  = 
I(u)I(y)I*(u  +  v ),  where  u  and  v  are  both  spatial  frequencies.  The  bispectra  of  the  object,  PSF,  and  speckle  image 
are  therefore  related  by  Bi(u,  v)  =  Bo(u,  v)Bh(u,  v).  It  can  be  shown  that  the  expected  bispectrum  (Bh(u,  v))  of  the 
atmospheric  PSF  has  zero  phase  (this  will  be  discussed  below),  so  the  ensemble  average  speckle  bispectrum  has  the  same 
phase  as  the  object  bispectrum.  From  the  bispectrum  phase,  the  object’s  Fourier  phase  may  be  recovered  up  to  an  overall 
linear  phase,  which  corresponds  to  a  translation  of  the  reconstructed  image. 


To  implement  this  approach,  a  sequence  of  speckle  images  is  collected  and  an  estimate  formed  of  the  bispectrum  of 
each  one.  In  the  presence  of  noise,  an  unbiased  estimator  of  the  bispectrum  is  given  by 

B,(u,  v )  =  I(u)I(v)I*(u  +  v)  -  \I(u)\2  -  \I(v)\2  -  I l(u  +  v)\2  +2K  +  3 npixa2,  (1) 

where  K  is  the  number  of  photoelectrons,  «p;x  is  the  number  of  pixels,  and  a2  is  the  read  noise  variance  in  each  pixel  in 
photoelectrons.  The  bispectra  of  the  individual  frames  are  averaged  together  to  form  an  ensemble  estimate. 

To  see  how  the  object  may  be  recovered,  consider  the  relationship  between  its  Fourier  phase  and  the  bispectrum,  given 
by  arg  Bo(u,v)  =  arg  O  (u)  +  arg  O  (v)  —  arg  O  (u  +  v) .  Thus  if  the  phase  is  known  at  u  and  v ,  it  may  be  recovered  at  u  +  v 
by 

arg  0(u  +  v)  =  arg  0(u)  +  arg  O(v)  —  ar gBo(u,  v).  (2) 

To  get  this  process  started,  observe  that  the  dc  term  is  always  real,  so  arg  0(0)  =  0.  The  choice  of  the  first  off-dc  Fourier 
components  governs  the  translation  of  the  reconstructed  image,  which  is  arbitrary.  From  this  point,  repeatedly  applying 
Eqn.  2  yields  the  entire  phase  map. 

Having  recovered  the  object’s  Fourier  phase,  one  needs  only  the  amplitude  1 0(u)\  to  obtain  a  complete  reconstruction. 
This  may  be  achieved  by  making  observations  of  a  star  under  similar  atmospheric  conditions  and  using  these  to  estimate 
the  speckle  transfer  function  (\H\2).  Then 

\0\  =  ((|/|2)/(|tf|2))1/2-  (3) 


2.2  Bispectrum  estimation  with  a  partitioned  aperture 

The  effectiveness  of  the  bispectrum  technique  hinges  on  the  convergence  of  the  bispectra  of  observed  atmospheric  PSFs  to 
the  purely  real  ensemble  average.  Taking  the  pupil  amplitude  to  be  unity  and  recalling  that  H(u)  is  the  autocorrelation  of 
the  pupil  complex  amplitude,  the  bispectrum  of  a  wavefront  with  phase  0  in  the  pupil  is 

Bh(u,v)=  /  exp(/(0(iu)  —  <p{w  +  u)))  dw 

Jp(0)np(u) 

x  /  exp(/(0(u;)  —  cp(w  +  v)))  dw  (4) 

Jp(0)np(v) 

x  /  exp +  u  +  w)  —  dw, 

Jp(0)np(u+v) 

where  P(z)  is  the  set  of  points  in  the  pupil  shifted  by  z.  As  the  number  of  atmospheric  cells  in  the  pupil  increases  (due  to 
increasing  D  or  decreasing  ro),  the  integral  product  yields  a  larger  number  of  uncancelled  phasor  contributions  to  Bh  (u  ,  v) . 
This  is  the  source  of  baseline  redundancy  noise. 

Suppose  that  instead  of  a  single  full  pupil,  speckle  images  are  collected  simultaneously  from  several  subregions  of  the 
aperture.  If  the  aperture  partition  is  constructed  in  such  a  way  that  the  number  of  points  in  each  region  having  a  given 
vector  separation  is  small,  then  the  redundancy  noise  in  the  bispectrum  estimate  will  be  reduced. 

Let  hj  (x)  be  the  PSF  associated  with  the  &th  subregion  and  let  ij  (v)  =  ( hj  *6>)(v)  be  the  corresponding  speckle  image. 
To  perform  a  bispectrum  recovery  from  this  information,  form  a  bispectrum  estimate  Bik(u ,  n)  for  each  region.  Then  in 
the  ensemble  average,  the  phase  of  each  region’s  bispectrum  will  match  the  object’s  bispectrum  phase: 

arg(B/jt(w,  v))  =  <xgB0(u,  v).  (5) 


Note  that  in  general,  each  region  will  have  a  different  diffraction  cutoff,  so  for  some  k  the  phase  on  the  left-hand  side  of 
Eqn.  5  could  be  undefined  for  some  u  and  v. 


A  simple  way  of  combining  the  bispectra  of  the  subregions  into  a  single  bispectrum  estimate 
pie  variance  of  the  estimates  and  use  that  to  form  a  weighted  mean  bispectrum.  Let  o2{u,  v) 
of  Bjk(u,  v).  Then  let 


^syn(^>  *0  —  ^  ' 


Bik(u,  V) 
crk(u,  v) 


/£ 


1 


ak(u,  v) 


is  to  calculate  the  sam- 
be  the  sample  variance 

(6) 


be  the  synthesized  bispectrum  estimate.  Note  that  this  is  weighting  by  standard  deviation,  not  by  variance  as  one  might 
expect.  In  simulation  results,  the  standard  deviation  was  found  to  produce  a  better  result. 
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Figure  1.  The  logarithm  (base  10)  of  the  bispectrum  SNR  in  the  subplane  u  =  (x,  0),  v  =  (0,  y).  On  the  left,  the  bispectrum 
is  formed  from  the  full  pupil.  On  the  right,  the  bispectrum  estimate  is  formed  from  an  aperture  divided  into  three  concentric 
annuli.  The  profile  was  scaled  to  yield  D/r$  =  16. 
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Figure  2.  The  bispectrum  SNR  for  the  full  and  partitioned  aperture  with  D/r o  =40. 


3.  SIMULATION  RESULTS 


3.1  Simulation  description 

To  test  the  performance  of  this  technique,  we  simulated  a  system  with  a  3.6  m  aperture  collecting  data  at  30  Hz  in  the  I  band. 
The  atmosphere  is  modeled  using  the  Maui-3  profile,  scaled  to  yield  various  turbulence  strengths.  The  atmosphere  is 
broken  into  four  slabs  to  capture  the  effects  of  slew  and  wind.  The  FPAs  are  256  x  256  CCDs  with  sky  foreground  and 
photon  noise  modeled,  but  no  read  noise.  Unless  otherwise  noted,  the  results  shown  are  for  an  aperture  partitioned  into 
three  concentric  annuli,  a  sidereal  slew  rate,  a  ground  wind  speed  of  4  m/s,  and  one  second’s  worth  of  data  (30  frames). 

3.2  SNR  results 

For  conventional  full-aperture  speckle  imaging,  various  authors  [3-5]  have  derived  analytic  approximations  of  the  SNR 
of  the  bispectrum  estimate.  These  calculations  are  based  on  the  assumption  that  the  aperture  contains  a  large  number  of 
statistically  independent  atmospheric  cells,  which  is  approximately  true  in  the  mid-frequency  region  with  a  full  aperture. 
However,  this  does  not  hold  for  the  annular  subapertures  that  are  of  interest  here.  Instead,  we  depend  on  simulations  to 
provide  measurements  of  the  bispectrum  SNR. 

Fig.  1  shows  the  bispectrum  SNR  in  the  plane  u  =  (x,  0),  v  =  (0,  y)  when  D/ro  =  16.  The  left  side  is  obtained  from 
the  full  aperture.  The  right  side  is  the  result  obtained  via  Eqn.  6  from  the  partitioned  aperture.  On  the  left  side,  the  usual 
features  are  visible  [3]:  very  high  SNR  near  the  origin,  high  SNR  along  the  \\u\\  ~  0  and  ||v||  ~  0  axes,  and  a  plateau  of 


Full  aperture  Partitioned  aperture 

Figure  3.  Reconstructed  images  from  a  full  aperture  (left)  and  partitioned  aperture  (right)  with  D/ro  =  64  and  sidereal 
slew  rate.  Both  data  sets  were  generated  from  the  same  simulated  wavefronts.  The  target  brightness  was  scaled  to  yield 
visual  magnitude  of  1 ,  which  gives  a  peak  SNR  of  820  on  the  full  aperture  FPA. 


much  lower  SNR  out  to  the  diffraction  limit.  The  partitioned  aperture  has  a  similar  structure,  but  the  regions  of  high  SNR 
extend  further  out  from  the  axes.  A  “ring”  structure  is  visible  in  the  SNR  map  due  to  the  annular  partitions  of  the  pupil. 

In  Fig.  2,  the  turbulence  strength  has  been  increased  to  give  D/ro  =  40.  The  color  scale  matches  that  of  Fig.  1,  so  the 
SNR  is  clearly  lower  everywhere.  As  before,  the  partitioned  aperture  yields  higher  SNR,  especially  in  the  region  off-axis 
but  close  to  the  origin.  The  importance  of  this  region  to  the  quality  of  the  object  phase  reconstruction  is  evident  from 
Eqn.  2.  If  Bo(u ,  v)  is  noisy  when  \\u\\  /  0,  ||u||  /  0,  then  the  estimate  of  0(u  +  v)  obtained  from  0(u)  and  0(v)  will 
also  be  noisy. 

3.3  Image  results 

To  get  an  idea  of  how  the  improved  SNR  translates  into  image  quality,  we  now  look  at  a  series  of  comparisons  between 
full-aperture  and  partitioned-aperture  reconstructed  images.  In  each  case,  the  same  wavefronts  are  used  to  generate  the 
simulated  images  from  which  the  reconstructions  are  formed. 

Fig.  3  illustrates  the  case  of  D/ro  =  64  for  a  bright  (visual  magnitude  mv  =  1)  target.  Recall  that  partitioning  the 
aperture  both  boosts  the  OTF  and  reduces  the  baseline  redundancy  noise.  Boosting  the  OTF  would  result  in  less  noise 
amplification  in  the  higher  spatial  frequencies,  but  in  this  case  the  SNR  is  so  high  that  there  is  virtually  no  sensor  noise  to 
amplify.  On  the  other  hand,  there  is  a  great  deal  of  baseline  redundancy  noise,  and  by  reducing  it  the  partitioning  method 
produces  a  much  better  result. 

In  Fig.  4,  the  target  brightness  is  varied  from  1st  mag  to  10th  mag.  In  each  case  the  partitioned  aperture  yields  the 
superior  result.  This  is  somewhat  counterintuitive,  since  splitting  the  aperture  to  form  separate  images  means  that  each 
partitioned-aperture  speckle  image  will  have  lower  SNR  than  the  full-aperture  speckle  images.  However,  under  the  condi¬ 
tions  simulated  here,  this  is  not  sufficient  to  outweigh  the  advantages  of  the  partitioning  technique.  Presumably  this  would 
no  longer  be  the  case  if  the  sensor  read  noise  were  sufficiently  high. 

As  the  turbulence  strength  increases  (Fig.  5),  the  degree  of  baseline  redundancy  increases,  causing  the  full-aperture 
images  to  suffer.  Breaking  up  the  pupil  suppresses  this  source  of  error,  greatly  improving  the  reconstructed  image  quality. 

Fig.  6  shows  the  effect  of  slew  rate  on  image  quality.  As  the  slew  rate  increases,  the  phase  in  the  pupil  decorrelates  more 
quickly.  This  improves  the  performance  of  the  full-pupil  results  because  redundancy  noise  can  be  effectively  averaged  out 
using  independent  realizations  of  the  atmosphere.  It  has  less  effect  on  the  partitioned-aperture  setup,  which  enjoys  a  much 
lower  level  of  redundancy  noise. 

Finally,  we  consider  the  effect  of  the  changing  the  number  of  annuli  into  which  the  pupil  is  partitioned.  As  shown 
in  Fig.  7,  if  there  are  too  few  annuli,  then  there  is  enough  redundancy  in  each  annulus  to  yield  perceivable  noise  in  the 
recovered  image.  However,  as  the  number  of  annuli  increases,  the  available  light  is  split  between  more  and  more  FPAs, 
reducing  the  SNR  on  each  one.  For  the  system  modeled  here,  the  best  performance  is  obtained  with  3  or  4  annuli. 
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Figure  4.  Full  aperture  vs.  partitioned  aperture  as  a  function  of  target  brightness.  The  peak  pixel  SNR  on  the  full  aperture 
FPA  varies  from  13  at  10th  mag  to  820  at  1st  mag. 
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Figure  5.  Full  aperture  vs.  partitioned  aperture  as  a  function  of  turbulence  strength.  The  increase  in  redundancy  noise 
causes  an  obvious  deterioration  in  the  quality  of  the  full-aperture  results. 
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Figure  6.  Full  aperture  vs.  partitioned  aperture  as  a  function  of  slew  rate  with  mv  =  4.  As  slew  rate  increases,  the  measured 
atmospheric  realizations  become  completely  decorrelated,  beating  down  redundancy  noise. 
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Figure  8.  Pupil  masks,  (a)  Aperture  mask  plate  installed  in  the  sensor,  (b)  Shack-Hartmann  images  of  each  subaperture. 
Note  that  these  were  made  with  an  internal  source,  so  the  secondary  obscuration  is  not  visible  in  the  smallest  mask. 


4.  EXPERIMENTAL  DESIGN 

We  conducted  an  experimental  test  of  aperture  partitioning  using  the  GEMINI  sensor  on  the  MSSS  1.6  m  telescope, 
which  is  capable  of  collecting  critically- sampled  imagery  at  over  200  Hz.  To  be  operationally  useful,  a  phase  mask  or 
mirror  arrangement  would  be  needed  to  partition  the  aperture  so  that  all  subapertures  could  be  measured  simultaneously. 
However,  as  an  initial  demonstration,  we  instead  used  a  series  of  pupil  masks  to  measure  each  subaperture  in  succession. 
This  is  a  valid  demonstration,  since  the  bispectrum  of  each  subaperture  is  estimated  independently.  However,  it  cannot  be 
applied  to  fast-moving  satellites,  since  the  object  pose  changes  noticeably  in  the  time  required  to  switch  masks. 

A  set  of  four  pupil  masks  were  made  in  a  steel  plate  (see  Fig.  8a)  and  mounted  on  a  motorized  linear  stage  at  a  pupil  in 
the  system.  Switching  from  one  mask  to  another  can  be  done  from  the  sensor  console  and  takes  about  two  or  three  seconds. 
The  masks  include  the  full  pupil  and  three  annular  subapertures.  The  outermost  mask  spans  the  entire  1.57  m  pupil,  the 
middle  annulus  has  an  outer  diameter  of  1.05  m,  and  the  inner  annulus  is  0.53  m  across.  Within  the  inner  annulus  is  the 
central  obscuration  of  the  telescope,  which  is  0.36  m  across.  Thus,  the  inner  annulus  has  a  thickness  of  about  8.5  cm  and 
the  middle  and  outer  annuli  are  about  26  cm  thick.  Wavefront  sensor  images  of  each  subaperture  are  shown  in  Fig.  8b. 

In  order  to  get  a  wide  range  of  conditions,  data  was  collected  over  two  days  (once  in  mid-morning  and  once  in  mid¬ 
afternoon)  and  three  nights  (all  evening  terminator).  Stars  and  binary  pairs  were  collected  at  elevation  angles  between  20 
and  85  degrees.  A  variety  of  filters  were  used,  typically  in  the  700-900  nm  range.  For  each  data  set,  seeing  was  estimated 
from  the  short-exposure  imagery  using  the  technique  of  von  der  Lfihe  [6] .  Median  daytime  seeing  was  3  cm  and  median 
nighttime  seeing  was  8  cm. 

Fig.  9  shows  data  collected  in  each  mask.  It  is  easy  to  see  the  reduced  resolution  in  the  data  from  the  middle  and, 
especially,  the  inner  subapertures  due  to  the  smaller  pupil.  Also,  the  data  from  the  inner  subaperture  clearly  has  much 
lower  SNR  because  the  collection  area  for  this  mask  is  much  smaller  than  for  the  others.  This  illustrates  the  point  that 
although  splitting  the  pupil  into  more  subapertures  would  reduce  redundancy  noise,  each  one  would  also  suffer  more  from 
sensor  noise  (read  noise  and  shot  noise)  due  to  the  reduced  signal  level. 


5.  EXPERIMENTAL  RESULTS 

To  compare  the  relative  performance  of  full-aperture  and  partitioned-aperture  speckle  imaging,  we  generated  quasi-exper- 
imental  extended-target  data  from  the  collected  stars  (although  binary  stars  were  collected,  they  were  not  analyzed  for  this 
study).  This  was  accomplished  by  convolving  each  star  frame  ( h )  with  a  TASAT  rendering  of  SEAS  AT  ( o ),  normalized 
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Figure  10.  Reconstructions  for  two  daytime  data  sets. 


to  °i  =  1-  One  can  see  that  this  does  not  have  valid  photon  statistics,  so  Gaussian  noise  with  variance  ( o  —  o 2)  *  h  was 
added  to  each  image  to  approximately  corrects  the  statistics. 

To  form  a  reconstruction  from  the  three- subaperture  data,  the  bispectrum  was  estimated  from  the  first  N  frames  of  each 
of  the  three  masks.  The  object  Fourier  phase  is  recovered  from  the  weighted  average  of  the  bispectra. 

To  form  a  reconstruction  from  the  full-pupil  data,  three  groups  of  N /3  consecutive  frames  were  used.  The  groups  were 
spaced  several  seconds  apart.  The  purpose  of  this  is  to  ensure  that  the  partitioned-aperture  method  does  not  derive  some 
advantage  from  the  fact  that  the  atmosphere  has  time  to  decorrelate  while  the  pupil  mask  in  the  sensor  is  being  changed. 

Some  example  reconstructions  are  shown  in  Figs.  10  and  11.  For  most  data  sets,  aperture  partitioning  produces  a 
visually  preferable  result.  In  order  to  better  quantify  this,  a  second  set  of  extended  data  was  created  using  a  resolution 
target  as  the  pristine  object.  Then  each  reconstruction  was  visually  assessed  to  determine  the  minimum  resolvable  feature 
size  in  pixels.  The  results  of  this  procedure  for  each  data  set  are  plotted  in  Fig.  12. 

Two  trends  are  apparent  in  these  results.  First,  severe  seeing  conditions  (smaller  ro)  consistently  favor  aperture  par¬ 
titioning.  In  particular,  aperture  partitioning  produced  a  better  result  for  all  but  one  of  the  daylight  data  sets.  This  is 
consistent  with  the  view  that  there  is  a  tradeoff  between  redundancy  noise  and  sensor  noise.  Using  a  partitioned  aperture 
reduces  redundancy  noise,  but  by  dividing  the  light  between  three  separate  focal  planes,  it  experiences  more  sensor  noise. 
As  ro  drops  (as  in  the  daytime,  for  example),  redundancy  noise  becomes  dominant.  But  in  good  seeing  conditions,  the 
lower  sensor  noise  in  the  full-aperture  data  produces  a  better  result. 

Second,  as  more  frames  are  processed  together,  the  two  methods  tend  to  approach  a  common  result.  This  is  because  as 
more  frames  are  used,  the  average  measured  bispectrum  approaches  the  true  ensemble  average  even  for  the  full-aperture 
case.  In  practice,  the  number  of  frames  that  can  be  processed  together  is  dictated  by  the  rate  at  which  the  object  (the 
satellite)  is  changing  pose. 
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Figure  11.  Reconstructions  for  two  nighttime  data  sets. 
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Figure  12.  Difference  in  log2  of  minimum  resolvable  feature  size  vs  ro  for  (a)  20  frames  and  (b)  320  frames.  Negative  values 
indicate  cases  in  which  aperture  partitioning  produced  the  better  result.  Its  advantage  tends  to  diminish  as  ro  increases. 


6.  CONCLUSIONS 


We  have  described  a  technique  for  improving  speckle  imaging  performance  by  partitioning  the  telescope  aperture.  In 
both  simulations  and  experiments,  we  find  that  when  D/ro  1  and  the  slew  rate  is  fairly  low,  it  greatly  outperforms 
conventional  full-aperture  imaging.  Unlike  aperture  masking  interferometry,  this  advantage  is  maintained  even  for  dim 
targets. 
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